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ABSTRACT 



Context. The interpretation of water line emission from existing observations and future HIFI/Herschel data requires a detailed 
knowledge of collisional rate coefficients. Among all relevant coUisional mechanisms, the rotational (de)excitation of HtO by H2 
molecules is the process of most interest in interstellar space. 

Aims. To determine rate coefficients for rotational de-excitation among the lowest 45 para and 45 ortho rotational levels of H2O 
colliding with both para and ortho-Ha in the temperature range 20-2000 K. 

Methods. Rate coefficients are calculated on a recent high-accuracy H2O-H2 potential energy surface using quasi-classical trajectory 
calculations. Trajectories are sampled by a canonical Monte-Carlo procedure. H2 molecules are assumed to be rotationally thermalized 
at the kinetic temperature. 

Results. By comparison with quantum calculations available for low lying levels, classical rates are found to be accurate within a 
factor of 1-3 for the dominant transitions, that is those with rates larger than a few 10^'^cm^^s"'. Large velocity gradient modelling 
shows that the new rates have a significant impact on emission line fluxes and that they should be adopted in any detailed population 
model of water in warm and hot environments. 

Key words, molecular data - molecular processes - ISM: molecules 



1. Introduction 



Since its discovery in interstellar space (Cheung et al. I1969I I. 
water vapour has been detected in a great variety of astronom- 
ical objects using both Earth-based and, foremost, spacecraft 
observations. The Infrared Space Observatory (ISO), in par- 
ticular, has revealed the ubiquity of water in the interstellar 
(ISM) and circumstellar (CSM) media (for a recent review see 
Cernicharo & Crovisier 2005). The more recent Submillimetre 
Wave Astronomy Satellite (SWAS) and Odin missions have also 
measured gaseous water in a wide variety of sources using 
the single ground-state rotational transitions of ortho-H20 at 
557 GHz (e.g. Melnick & Bergin |2005l Hjalmarson et al. |2b03l l. 
Finally, water is also detected at radio and (sub)millimeter wave- 
lengths through maser transitions which are commonly associ- 
ated with star-forming regions (e.g. Cernicharo et al. |1990|[1999l 
Goddi et al. 1 20071 and references therein or extraglactic sources 
(e.g. Cernicharo et aI. I2006l Kondratko et aI. I2006l and references 
therein). 

The abundance of water varies largely in cold and warm me- 
dia. In the former, because of the freezing of water onto the dust 
grains, water does not exceed ~ 10 (Bergin & Snell 120021 
Poelman et al. 120071 and references therein). But in warm re- 
gions, water can become the most abundant molecule after H2, 
with an abundance of ~ 10 In star forming regions, this oc- 
curs both because of the release in the gas phase of the grain ices 
and because of endothermic reactions efficiently forming water 
in the gas phase (e.g. Ceccarelli et al. il996t Doty & Neufeld 
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II 9971 and references therein). Under these conditions, water be- 
comes a crucial molecule in the thermal and chemical balance. 
In fact, owing to its large dipole, water lines will dominate the 
cooling (and sometime even heating) of the gas in a large range 
of gas densities and temperatures (e.g. Neufeld & Melnick[1991j 
Ceccarelli et al. |1996t Cernicharo et al. 120061 and references 
therein). Also, since oxygen will be mostly locked in water 
molecules at temperatures larger than about 200 K, the chem- 
ical abundance of more complex and less abundant molecules 
will depend on the water abundance. 

For all these reasons, water is a key molecule in space, and 
its observation and the determination of its abundance are criti- 
cal for many theoretical aspects. Indeed, observing and measur- 
ing the water lines and their abundance is an important task for 
the high resolution spectrometer Heterodyne Instrument for the 
Far Infrared (HIFI) on board the ESA funded Herschel satel- 
lite, which will be launched in 2008. However, as stated by 
Cernicharo & Crovisier ( 120051 ). "little new information will be 
obtained if collisional rates adapted to the temperatures of the 
clouds in the ISM and CSM are not available". 

Collisional rate coefficients are indeed essential in the de- 
scription of the energy exchange processes responsible for 
molecular line formation in astronomical environments. Spectral 
features such as maser emission are produced in low-density 
conditions far from thermodynamic equilibrium and through 
a complex competition between radiative and collisional pro- 
cesses, including infrared emission from dust (e.g. Poelman et 
al. l2007l l. In the ISM and CSM, the main colliding partners are 
hydrogen molecules, hydrogen and helium atoms, and electrons. 
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In general, however, H2 is the dominant form of hydrogen and 
as it is more abundant than He (by a factor of ~ 5) and elec- 
trons (by ~ 4-8 orders of magnitude), H2 is usually the dom- 
inant exciting species. At typical ISM and CSM temperatures 
(T < 2000 K), the dominant energy exchange process is thus 
rotational (de)excitation by H2, although vibrational excitation 
cannot be neglected in the highest temperature and density re- 
gions. 

To the best of our knowledge no experimental state-to-state 
data is available for H2O-H2 except for the total relaxation rate 
of the water bending mode (see Faure et al. 2005a and references 
therein). On the theoretical side, quantum close-coupling (CC) 
state-to-state calculations for the rotational excitation of H2O by 
H2 and He have been performed for kinetic temperatures and 
energy levels lower than 140 K and 2000 K, respectively (Green 
et al. 1993; Phillips et al. [T996l Dubernetet al. 2002; Grosjeanet 
al.|2003; Dubernet et al. l2006l l. Hence, current models of water 
emission in warm and hot environments rely exclusively on rates 
for excitation by He atoms, usually scaled by the reduced mass 
ratio (jUH,o-Ho/A'HoO-He)''^ (e.g. Poelman & Spaans'2005). It has 
been shown by Phillips et al. ( 1996), however, that excitation by 
para-H2 in its 7=0 level is not too different from excitation by 
He atoms (with most rates differing by a multiplicative factor 
of 1-3) but that excitation by H2 in excited rotational levels is 
significantly different. This mainly reflects the importance of the 
long-range interaction between the H2 quadrupole moment and 
the dipole of H2O (PhilHps et al. 119961 ). Rates for excitation of 
H2O by H2 at high temperatures and/or for high rotational levels 
are therefore urgently needed for the interpretation of existing 
and future water emission spectra. 

In principle, scattering calculations based on the quantum 
CC method can provide an absolute accuracy of a few percent 
for a given potential energy surface (PES). As an illustrative 
example, the recent study of Gilijamse et al. (12006 !) has shown 
an excellent agreement between experimentally measured cross 
sections and CC calculations based on ab initio PES. The ma- 
jor drawback of the quantum CC method is its computational 
cost which increases dramatically with the number of coupled 
channels, i.e. with the total energy and the number of degrees 
of freedom. Very recently, Dubernet et al. (120061) provided new 
rotational rate coefficients for H2O-H2 at temperatures below 
20 K using fully converged CC calculations based on the recent 
high-accuracy PES of Faure et al. (12005a.) . The main objective 
of their study was to test the influence of the new PES on the 
scattering calculations and to make comparisons with previous 
results based on the PES of Phillips et al. ( 1994). A significant 
re-evaluation of rates was observed, especially for para-H2, with 
differences up to a factor of 3 at 20 K. Full CC calculations, 
augmented with coupled-states and infinite order sudden (lOS) 
calculations, are currently in progress at higher collision ener- 
gies with the objective of obtaining quantum CC rates for kinetic 
temperatures and rotational levels up to ~ 2000 K and making 
comparisons with pressure broadening experiments (Dubernet 
et al., in preparation). 

The present study is motivated by the possibility of per- 
forming high temperature quasiclassical trajectory (QCT) cal- 
culations as an alternative to quantum mechanical calculations. 
In contrast to CC calculations, the computational time required 
for QCT calculations decreases as the collisional energy in- 
creases and such calculations are thus particularly adapted to 
high temperatures where purely quantum effects are usually neg- 
ligible. QCT calculations have been shown to accurately repro- 
duce quantum results provided that the relevant transition prob- 
abilities are large enough. In particular, Faure et al. ( 120061 ) re- 



cently observed a good agreement between classical and quan- 
tum rates for the rotational de-excitation of H2O by H2 at a sin- 
gle temperature of 100 K. The average accuracy of classical re- 
sults was found to be better than a factor of 2 for the transitions 
with the largest probabilities, i.e. those with a rate greater than 
~ 10""cm^s"'. The classical approach is thus expected to pro- 
vide a good compromise between accuracy and computational 
effort at temperatures where converged quantum CC calculations 
are extremely expensive (typically T >300 K). The next section 
describes details of the QCT calculations. Rate coefficients are 
presented in Sect. 3. A first application of these rates is given in 
Sect. 4. Conclusions are drawn in Sect. 5. 

2. Quasi-classical trajectory method 

2.1. Potential energy surface and general procedure 

All calculations presented below were performed with rigid 
molecules using the vibrationally averaged PES of Faure et al. 
( I2005al l. where full details can be found. The angular expansion 
of the PES was obtained as a least square fit of thousands of ab 
initio values over 149 angular functions which are specifically 
adapted to quantum calculations and provide analytical deriva- 
tives for the classical equations of motion. In order to reduce 
computational time in integrating these equations, a subset of 
only 45 angular functions was retained. This selection was found 
to reproduce the ab initio values within a few cm in the Van 
der Waals minimum region. The inclusion of additional angu- 
lar functions was tested through trajectory analysis and found to 
be negligible. The water and hydrogen geometries were taken 
at their effective rotational values, that is those corresponding to 
the spectroscopically determined rotational constants. 

As inelastic (nonreactive) collisions cannot interconvert the 
ortho- and para-forms of either species, QCT calculations were 
done separately for the four nuclear spin combinations, namely 
para-H20-para-H2, para-H20-ortho-H2, ortho-H20-para-H2 
and ortho-H20-ortho-H2. Classical trajectories were run using 
the Monte-Carlo procedure described in Faure et al. (^2006). 
The rate coefficients were computed by a canonical sampling 
of the initial Maxwell-Boltzmann distribution of collisional en- 
ergy at four kinetic temperatures, namely T =100, 300, 1000 
and 2000 K. Some test calculations were also performed at 50 K 
(see Fig. 1). The classical equations of motion were numerically 
solved using an extrapolation Bulirsch-Stoer algorithm. The to- 
tal energy was conserved up to six digits, i.e. within ~0.1 cm"'. 
The maximum impact parameter /Jmax, defined as the small- 
est impact parameter for which 1000 trajectories produced no 
change in the rotational state of H2O, was found to range be- 
tween 10 and 16 a.u. at the selected temperatures. Batches of 
10,000 trajectories were run for each of the lowest 44 para and 
44 ortho excited rotational levels of H2O, resulting in a total 
of 7,040,000 trajectories. These calculations required approxi- 
mately 150 hours of computer time on AMD Opteron quadri- 
processors for each temperature and for each of the four nuclear 
spin cases. 

2.2. Initial state selection 

For each molecule, the classical rotational angular momentum 
was taken as [J(J + 1)]'^^, where J is the rotational quantum 
(integer) number, using the correspondence principle. For wa- 
ter, initial values of the classical angular momentum components 
were simply assigned to the pseudoquantum numbers Kg and 
(projection of J along the axis of the least and greatest moments 
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of inertia, respectively). This trivial initial selection was found 
to reproduce exact quantum eigenvalues within a few cm"' for 
the lowest levels (see Table 1 in Faure et al. 2006) and within 
a few percent for the highest. For H2, the initial para and or- 
tho rotational levels were weighted separately according to the 
Boltzmann distribution: 



Eh 



(272 + l)exp(-^) 
2:y,(272 + l)exp(-|^) 



(1) 



where J2 is the H2 angular momentum, k is the Boltzmann con- 
stant and the sum in the denominator extends only over even 
(odd) numbers for para (ortho) H2. We thus assumed that each 
nuclear spin species of H2 has a thermal distribution of rotational 
levels. This is obviously not necessarily the case in astronomical 
environments, in particular for high-72 levels. However, from the 
H2 Einstein coefficients tabulated by Wolniewicz et al. (I1998I I 
and the H2-He rates reported by Flower et al. ( 1998) , we have 
estimated the He critical densities to be below ~ 10^ cm up to 
J2 - 5, that is up to H2 levels lying more than ~2000 K above the 
ground states. These levels are therefore likely to be at or close 
to local thermodynamic equilibrium (LTE) at typical ISM and 
CSM densities. Observationally, the detections of H2 in shocks 
by Lefloch et al. ( 120031) and Neufeld et al. (I2006I I also suggest 
that the H2 rotational populations are close to LTE for J2 as high 
as 9. As the dependence of the H2O rates on the H2 rotational 
level is modest (see Sect. 3), the assumption of rotationally ther- 
malized H2 molecules seems reasonable. Furthermore, this as- 
sumption has the advantage that the rate coefficients rigorously 
satisfy the principle of detailed balance (see below). 

2.3. Final state selection and rate calculation 

In the case of an asymmetric top molecule, the standard bin his- 
togram method used to assign the final classical action has been 
shown to be ambiguous (Faure & Wiesenfeld i2004) . However, in 
the particular case of water-like molecules, i.e. asymmetric ro- 
tors with /) dipole along the axis of intermediate moment of iner- 
tia and //) C2V symmetry, the ambiguities can be solved by taking 
advantage of the fact that ortho/para transitions are not permitted 
(see Faure et al. l2006l l. For H2, no final assignment was required 
as rate coefficients must be summed over all possible H2 rota- 
tional levels. The rate coefficients reported in the present study 
are therefore defined, for each of the four spin combinations, as: 



R(ja ^ 4,) = ^p(/2)-R(^,/2 ^ /;-,/2)' 



(2) 



J2,A 



where a - Ka, Kc. According to the canonical Monte-Carlo pro- 
cedure described above, these rates were computed as: 



1/2 



8kT 



N_ 

'Mot 



(3) 



where ji is the reduced mass, N is the number of trajectories 
leading to a specified final level of water J'^, and Mot is the total 
number of 'physical' trajectoriefl typically between 9,000 and 
10,000. The Monte-Carlo standard deviation is: 



w,. ■/;,) _/ Mot-M 



1/2 



(4) 



' A small fraction (less than 10%) of trajectories are indeed rejected 
due to either non-conservation of energy (failure of numerical integra- 
tion) or unphysical final values of K„,K^. 



In the following, statistical errors are defined as two standard 
deviations which, in practice, correspond to a typical statistical 
accuracy of 10-30%. 

2.4. Detailed balance 

The principle of detailed balance (microsopic reversibility), 
which is the consequence of the invariance of the interaction un- 
der time reversal, states that the rate coefficients, as defined in 
Eq. (2), must obey the relation: 



R{J'^, ^ Ja, T) 



(2J+1) 
(27' + 1) 



exp 



Ej' -E, 



kT 



R(Ja^J'^,,T).(5) 



As shown by Faure et al. (120061 ), however, QCT calculations 
significantly overestimate the excitation (endoergic) rates in the 
case of H2O transitions whose energies are larger or comparable 
to the kinetic temperature. This 'threshold effect' results from 
the rotational degrees of freedom of H2O that are classically ac- 
tive at energies below quantum thresholds. The competing effect 
induced by the rotational degrees of freedom of H2 was also ob- 
served, to a lesser extent, in the case of H2O transitions whose 
energies are lower than the kinetic temperature. Such effects 
were also reported for other hydrogenic systems (see Mandy 
& Pogrebnya 2004 and references therein). As a result, detailed 
balance is not rigorously satisfied at the classical level and ex- 
citation rates can be ovestimated by an order of magnitude. The 
recommended practical procedure is therefore to employ QCT 
calculations for de-excitation (exoergic) transitions and to apply 
Eq. (5) to obtain the reverse, endoergic, rates. We note, however, 
that detailed balance was actually satisfied within error bars for 
a number of low energy transitions at the highest temperature^. 

2.5. Advantages, shortcomings and accuracy of the classical 
method 

The shortcomings of the classical approach are well known and 
are usually marked at low kinetic temperatures where purely 
quantum effects are significant. In particular, between ~ 1 and 
100 K where the Van der Waals interactions are dominant, quan- 
tum rates exhibit oscillatory temperature dependence due to 
the presence of quantum resonances (of shape and/or Feshbach 
type). In this regime, classical calculations are obviously ques- 
tionable. Amplitudes of these oscillations are however usually 
modest and not much larger than the accuracy of QCT calcula- 
tions. We note, for instance, that the recent study of Wernli et al. 
(120071) showed a quite good agreement between QCT and quan- 
tum rates down to 10 K. A greater limitation of the classical ap- 
proach is the quantization of the rotational degrees of freedom, 
as explained above. Finally, a major disadvantage of classical 
mechanics is its inability to reproduce small probabilities, the 
so-called classically forbidden transitions corresponding to rates 
lower than ~10"'^ - 10 '^cm^s^'. Owing to all these intrinsic 
limitations that add to the Monte-Carlo statistical deviations, an 
accuracy of greater than 10% is difficult to achieve by QCT rate 
calculations. The recent work of Faure et al. (2006) has shown, 
however, that classical rates for H2O-H2 are accurate within a 
multiplicative factor of 1-3 for low lying levels. Precision of 
-30% was even observed for the greatest rates. Such accuracies 



- It should be noted that 'effective' rates which depend on the initial 
rotational levels of H2 and which are summed over the final H2 rota- 
tional levels (as defined in Eq. (1) of Dubernet et al. 12006 ) do not obey 
detailed balance since they do not involve a complete thermal average. 
These effective rates were reported in the previous quantum studies. 
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are still much lower than the few percent precision of quantum 
CC calculations. 

However, a major advantage of classical trajectories is to al- 
low the computation of de-excitation rates for high initial values 
of the rotational angular momenta. In particular, the calculation 
of rate coefficients as defined in Eq. (2) requires the inclusion of 
all significantly populated levels of H2. At 2000 K and assuming 
LTE, the most populated level is ^2=4 (3) for para(ortho) H2 but 
levels up to 72=14 (13) were in fact included in our QCT calcu- 
lations. At the quantum CC level, a calculation with ^2 > 4 is 
cuiTently not feasible. Thus, classical mechanics proves particu- 
larly suited to collisions between 'hot' molecules where standard 
quantum methods are not applicable. 



3. Quasi-classical trajectory results 

3.1. Comparison to quantum data 

In Fig. 1 , four sets of rate coefficients for rotational de-excitation 
of the ortho-H20 3o3 level by para and ortho-H2 are plotted as 
a function of temperature: i) current QCT data between 50 and 
2000 K, ii) quantum data of Phillips et al. d 19961 1 be tween 20 and 
140 K, Hi) H20-He quantum data of Green et al. ('1993') multi- 
plied by the reduced mass ratio (a/h,o-h. //"H^o-He)'^^ - 1.348 
between 20 and 2000 K, iv) quantum data of Dubernet et al. 
(I2OO61) at 20 K augmented with unpublished data for para-H2 be- 
tween 20 and 1500 K. These latter unpublished rates have been 
computed using Eq. (2) with the sum extending over J2 = 0, 2, 4 
and with a scaling procedure for estimating the rates for exci- 
tation by H2(72 - 4) from those by H2(72 - 2fl We have 
selected the 3o3 level because it is the highest initial state in- 
cluded by Phillips et al. (119961 ). Firstly, it can be observed that 
there is an overall good agreement (within a factor of 1 -3) be- 
tween the present rates and those of both Phillips et al. d 19961 ) 
and Dubernet et al. while large differences (up to 2 orders of 
magnitude) are found with the scaled H20-He rates, especially 
so for ortho-H2 at T <100 K. It should be noted in particular 
that at low temperatures, the transition 3o3 - lio has a larger 
rate than 3o3 - 221 for para-H2 while the reverse holds for ortho- 
H2. Second, in the case of para-H2, the agreement between the 
present classical calculations and the quantum calculations of 
Dubernet et al. is excellent above 100 K, where quantum data are 
within the eiTor bars of the QCT data. At and below 100 K, the 
agreement is not as good, reflecting the limitations of the classi- 
cal approach (see Sect. 2.5), but differences do not exceed a fac- 
tor of 2. Third, the influence of the new H2O-H2 PES is modest, 
except for para-H2 at low temperature. Classical rates for ortho- 
H2 are thus found to be in very good agreement with those of 
Phillips et al. (|1996 ). This small impact of the PES on the ortho- 
H2 rates was previously observed by Dubernet et al. ( 120061 ). 
Fourth, we observe that the differences between the present clas- 
sical rates and the scaled H20-He rates decrease with increasing 
temperature. This reflects the fact that at high collision energy, 
the scattering process becomes dominated by kinematics rather 
than specific features of the PES. Finally, above 300 K, the dif- 
ference between para and ortho-H2 rates becomes minor. This 
clearly shows that the dependence of the rates on the H2 rota- 
tional level is weak, provided that ^2 ^ 1- We have actually 
checked that differences between de-excitation rates of the 3o3 
level by H2(72 - 1) and H2(72 = 6) do not differ by more than 
four Monte-Carlo standard deviations, i.e. by less than a factor of 



2. This finding is consistent with the quantum results of Phillips 
et al. ( 119961 ) who found similar rates for excitation of para-H20 
byH2(72 = l)andH2(72 =2). 
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Fig. 1. Rate coefficients, Eq. (2), for rotational de-excitation of 
the ortho-H20 state 3o3 as functions of temperature for para- 
H2 (upper panel) and ortho-H2 (lower panel). Squares denote 
quantum results by Phillips et al. (1996 ); stars refer to the quan- 
tum calculations of Dubernet et al.; dashed lines denote quan- 
tum (scaled) results of Green et al. (,1993i ) for H20-He; cir- 
cles with error bars (95% statistical confidence level) give the 
present QCT results. See text for details. Note that the transition 
3o3 — » loi is not plotted for clarity. 



3.2. Accuracy and propensity rules 

The above observations obviously hold for other initial levels 
and, from numerous comparisons, we estimate the QCT results 
to be accurate within a factor of 1-3 for the largest rates, i.e. 
those larger than a few 10^'^cm^s"^ as observed in Faure et al. 
(120061 ). For smallest rates, errors larger than a factor of 3 can be 
encountered, illustrating the inability of classical mechanics to 
reproduce small probability events. As a result, rates less than ~ 
lO^^^cm-'s which correspond to very poor statistics (zero or 
less than ~10 events over 10,000), were set to zero. This means 
that a significant number of transitions, the so-called classically 
forbidden or rare transitions, have null classical rates. Moreover, 
some transitions were found to have null rates at 100 K but 
significant rates at higher temperatures. Our strategy to tackle 
this problem was to retain transitions obeying the following two 
criteria: the rate must be /) larger than 10^'^cm^s"' at 300 K 
and above and //) larger than 10""cm^s"' at 2000 K. If one 
of the two criteria is not fulfilled, then the transition is clas- 
sified as a rare event and the corresponding (scaled) H20-He 
rate is adopted. It should be noted that these rare events corre- 
spond in most cases to transitions with A7 > 3 which are not 
expected to play a dominant role in the radiative transfer equa- 
tions (see Sect. 4). Moreover, despite the fact that no rigorous 
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selection rules hold for molecule-molecule collisions, we actu- 
ally observed the following propensity rules: 



A7 = 0,+l; AKa=Q,+l\ AKr^Q,+l, 



(6) 



which are consistent with the classical picture of a minimum re- 
orientation of the angular momentum of a strongly asymmetric- 
top moleculfl These rules are illustrated in Fig. 1 where the 
transition 3o3 - 2i2 is indeed found to be strongly favoured. 



3.3. Fitting procedure 

For use in astrophysical modelling, the transition rates RiJ',^, — » 
Ja, T) must be evaluated on a sufficiently fine temperature grid. 
For this, the rates were first least-square fitted over 100-2000 K 
by the analytic form used by Faure et al. ( 120041 : 



iogio^(/;- = 



(7) 



n=0 



where x = 7^-1/6 ^ typical fitting accuracy of a few percent 
was obtained. In the case of rare events with large error bars, the 
fitting method was also found to smooth the temperature depen- 
dence of the rates and thus remove iiTegularities caused by poor 
statistics. Note that T =50 K was investigated only for a few lev- 
els and was not retained in the fitting procedure. We then eval- 
uated the fitted rates at the following grid temperature values: 
100, 200, 400, 800, 1200, 1600 and 2000 K. From the H2O-H2 
data of Phillips et al. (1996), we observed the temperature de- 
pendence of de-excitation rates below 100 K to be very weak: 
de-excitation rates at 20 K are always within a factor of 2 of 
those at 100 K. This is clearly illustrated in Fig. 1. Below 100 K, 
a flat temperature dependence of the de-excitation rate was thus 
assumed, except for the lowest five para- and five ortho-HaO lev- 
els for which the values at 20 K were taken from the quantum 
CC results of Dubernet et al. (I2006I I. Such a flat low tempera- 
ture dependence actually reflects the influence of the deep po- 
tential well of the H2O-H2 PES (~ 235 cm ') which supports 
large shape resonances, in contrast to H20-He (~ 35 cm 
Obviously, the above low temperature extrapolation procedure 
introduces errors but these are believed not to exceed a factor of 
2-3 for the largest rates. The complete set of de-excitation rates 
among all 45 para- and ortho-H20 levels between 20 and 2000 K 
are provided as online material and will be made available in 
the BASECOL database (http://www.obspm.fr/basecol/). 
Excitation rates can be obtained from the detailed balance rela- 
tion, Eq. (5). 

In order to assess quantitatively the impact of the new 
H2O-H2 collisional rates on water emission, detailed radiative 
transfer studies adapted to various astronomical environments 
are required. In Sect. 4 below, this impact is investigated for a 
range of density and temperature conditions by neglecting all 
excitation mechanisms other than collisions. A first-order indi- 
cator is also provided by the impact of rates on critical densities. 
For a multi-level system, the critical density of a given partner 
(here para- or ortho-H2) is usually defined, in the optically thin 
case, given that the density at which the sum of the collisional 
de-excitation rates of a given level is equal to the sum of the 
spontaneous radiative de-excitation rates: 



ricriJa, T) 



(8) 



Table 1. Critical densities, in cm"^, for (scaled) He, para-H2 and 
ortho-H2, as functions of temperature and for a set of representa- 
tive H2O levels. Powers of 10 are given in parentheses. See text 
for details. 
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In the case of a slightly asymmetric-top such as H2CO, K„ and 
are further constrained by the value of J . 



When the colliding partner has a density n » ricr, collisions 
maintain rotational levels in LTE at the kinetic temperature while 
for densities n < n^r, deviations from LTE including population 
inversions are expected. Eq. (8) was computed with the present 
collisional rates for para-H2 and ortho-H2 and with the scaled 
H20-He rates of Green et al. (II9931 ). Einstein coefficients were 
taken from the Jet Propulsion Laboratory catalogue. Results are 
presented in Table 1 for 6 representative water levels, namely 
the upper level of the 557 GHz transition (1 10 - loi) and the up- 
per levels of maser transitions known to be collisionally pumped 
(Yates et al. [T997] ): the 183 GHz fine (3i3 - 220), the 380 GHz 
line (4i4 - 32i), the 325 GHz line (5 15 - 422), the 22 GHz fine 
(616 - 523) and the 321 GHz line (IO29 - 936). Note that maser 
transitions are expected to be quite sensitive to relatively small 
differences in collisional rates. As illustrated in Table 1, critical 
densities for H2 are decreased with respect to the scaled He val- 
ues by a typical factor of 2-5 at 100 K and by a factor of 1-3 at 
higher temperatures. These differences are still modest (with the 
exception of IO29) because collisional rates are summed over all 
possible downward transitions, thus averaging the impact of the 
new rates. Larger differences in the emission line fluxes are how- 
ever expected as individual state-to-state collisional rates, in par- 
ticular with ortho-H2, can exceed the scaled He values by more 
than an order of magnitude. This is demonstrated in Sect. 4. We 
conclude that using the new H2O-H2 collisional rates, non-LTE 
effects including population inversions will be quenched at lower 
H2 densities, especially in regions where the H2 ortho/para ratio 
is large. 



4. Modelling of water emission 

In this section we explore, by means of a LVG code, the im- 
pact of the new computed collisional coefficients on the theo- 
retical predictions of water line emission. The code is adapted 
from the one described in Ceccarelli et al. (120021 ). It refers to 
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a semi-infinite isodense and isothermal slab (plane-parallel ge- 
ometry). The code solves self-consistently the statistical level 
populations and the radiative transfer equations under the ap- 
proximation of the escape probability, assuming a gradient in the 
velocity field (Capriotti 1965 ). We consider the first 45 levels of 
both ortho- and para-H20 separately. The computed line spec- 
trum depends on a few basic parameters: the density of the col- 
liders (H2), the temperature of the emitting gas, the water column 
density and the maximum gas velocity. To have a relatively ex- 
haustive study, we explored a large parameter space. We varied 
the H2 density between 10"^ and 10'" cm"^; the kinetic tempera- 
ture between 20 and 2000 K; the water column density between 
10'" and lO'^cm"^, keeping the velocity equal to 1 km/s. Note 
that the velocity only enters in the line opacity, coupled with the 
H2O column density. In the present study we considered specif- 
ically the case of optically thin lines (the low H2O column den- 
sity case), and the case where the lines are optically thick (the 
high H2O column density case). In both cases, depending on the 
gas density, the collisional coefficients may play a major role in 
the predicted flux. Indeed, even in the case of strongly optically 
thick lines, the lines could be "effectively optically thin" if the 
levels are very sub-thermally populated, as is often the case for 
water lines. Finally, since the scope of this study is to understand 
the impact of the newly computed collisional coefficients in the 
line predictions, compared to the old ones, we did not consider 
the pumping of the levels by infrared and submillimetre radia- 
tion, which would unnecessarily complicate the problem. 

We compared the results of the LVG computations with three 
sets of collisional data: the scaled H20-He rates of Green et al. 
( fT993] l betw een 20 and 2000 K, the H2O-H2 rates of Phillips 
et al. d 19961 1 between 20 and 140 K and the present H2O-H2 
classical rates beween 20 and 2000 K. The results of this com- 
parison are reported in Fig. |2] where we show the line flux ra- 
tios for ortho-H20 at representative densities and temperatures. 
Collision rates with ortho-H2 only are considered to simplify the 
interpretation. The H2O column density is fixed at 10'^ cm"^. It 
can be observed that LVG fluxes based on the classical rates are 
increased up to a factor of 10 with respect to LVG fluxes based 
on the scaled H20-He rates of Green et al. (119931 ). Even larger 
differences were observed at lower temperatures. Second, LVG 
fluxes based on the classical rates are very similar to those based 
on the rates of Phillips et al. d 19961 1, with flux ratios close to 1 . 
These findings clearly reflect the differences in collision rates, 
such as those illustrated in Fig. 1. They also show that differ- 
ences in rates are not amplified within the radiative transfer equa- 
tions. In particular, it is worth noting that the typical increase of 
fluxes at 1000 K is less than a factor of 2 at the investigated 
densities. Similar results were obtained with para-H2, especially 
above 300 K where differences between para- and ortho-H2 rates 
are minor We note, however, that below 100 K, the classical 
para-H2 rates are not necessarily more accurate than the scaled 
H20-He rates. Third, we observe that the flux ratios rise steeply 
with increasing upper energies. This is expected since the critical 
densities are larger for higher levels and, therefore, the impact of 
collision rates, i.e. non-LTE effects, becomes more pronounced. 
Finally, when the temperature increases, the impact of collision 
rates is reduced, partly because the low lying levels become ther- 
malized and partly because the differences in rates decrease. It 
should be noted, however, that at densities higher than lO^cm"^, 
large flux ratios (~ 10) were still observed at 1000 K for high 
lying levels. 

We also tested the influence of the scaled H20-He rates 
employed as substitutes for H2 rates in the case of classically- 
forbidden transitions. These rates, which generally correspond 



to A7 > 3 and are lower than 10 cm^s were found to play 
a minor role, especially at low temperatures (T <300 K). This 
suggests that the radiative transfer equations are above all sen- 
sitive to the largest rates. These findings however merit further 
investigation. 

5. Conclusions 

We have computed classical rate coefficients for rotational de- 
excitation among the lowest 45 para and 45 ortho rotational lev- 
els of H2O colliding with both para and ortho-H2 in the temper- 
ature range 20-2000 K. Trajectories were run on a recent high- 
accuracy H2O-H2 PES using a canonical Monte-Carlo sam- 
pling. H2 molecules were assumed to be rotationally thermalized 
at the kinetic temperature. By comparison with quantum calcula- 
tions available for low lying levels, classical rates were found to 
be accurate within a factor of 1 -3 for the dominant transitions, 
i.e. those with rates larger than a few 10^'^cm^s"'. Large differ- 
ences with the scaled H20-He rates of Green et al. ( 1993) were 
observed and emission line fluxes were shown to be increased by 
up to a factor of 10. The present classical rates should therefore 
be adopted in any detailed population model of water in warm 
and hot environments. 

It should be noted, however, that quantum CC calculations 
are currently in progress (Dubernet and co-workers) and quan- 
tum rates will ultimately replace the classical rates. It will thus 
be possible in future studies to further test the sensitivity of wa- 
ter emission to uncertainties in the collision rates and to identify, 
for example, physical regimes where classical rates are of suffi- 
cient accuracy. We also stress that the QCT method is currently 
the only practical approach to investigate very high rotational an- 
gular momenta. For example, water in supergiant stars has been 
observed through levels with rotational quantum number / up to 
23 (e.g. Jennings & Sada J_998). In order to assess non-LTE ef- 
fects in such environments, classical calculations might provide 
precious information. Finally, we note that the first excited vibra- 
tional level of water lies only 2300 K above the ground state and 
it has been detected in several sources (e.g. Menten et al. 120061 ). 
Rate coefficients for the total vibrational relaxation of H2O by 
H2 molecules have been computed by Faure et al. (I2005al ) and 
some "7-resolved" cross sections can be found in Faure et al. 
(2005b ). Rovibrational state-to-state rates are however unknown 
and will be considered in future works. 
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Fig. 2. Emission line flux ratios from LVG modelling as functions of the upper energy levels of ortho-H20 for representative 
densities and temperatures. Three differents sets of collision rates were employed (see text for details). The H2O column density is 
fixed at 10'^ cm"^. Squares denote the ratios between LVG fluxes based on the classical rates and LVG fluxes based on the H2O-H2 
rates of Phillips et al. ( 1996 ). Diamonds and stars give the ratio between LVG fluxes based on the classical rates and LVG fluxes 
based on the scaled H20-He rate of Green et al. ( 119931 ). Diamonds correspond to transitions whose fluxes are larger than 1% of the 
total flux; stars coiTespond to the radio transitions at 22 and 380 GHz. 



